T 


LEVELS 


MO  16%1.1-A-e’s 


MECHANICAL  ENGINEERING  DEPARTMENT  ^ 


S.  SCHOOL  OF  ENGINEERING 

*4  THE  UNIVERSITY  OF  CONNECTICUT 


STORRS,  CONNECTICUT 

81  2  17  103 

w-4.  „ 


mmmmm  ■ 


i 

I 

t 

i 

i. 

EFFECT  OF  ADDITIVES  ON  THE 
LEAN  DETONATION  LIMIT  OF  KEROSENE  SPRAYS 
IN  AIR 

Final  Report 
by 

E.  K«  Dabora 


December  19 8C 


U.  S.  Army  Research  Office 
Grant  No. :  DAAG  29-78 -G0074 


THE  UNIVERSITY  OF  CONNECTICUT 
Department  of  Mechanical  Engineering 
Storrs ,  Connecticut  06268 


DISTRIBUTION 

Approved  {or  public  r«I*oM; 
Dfcttllwtlon  Unlimited 


mmmman 


ttrant. 


•m u 


FOREWORD 


This  report  summarizes  the  results  on  the  research  work  performed 
under  Army  Research  Office  Grant  No.  DAAG  29-78-G-0074  which  started  in 
March  1978  and  terminated  in  October  1979.  During  that  period  the 
University  of  Connecticut  Research  Foundation  provided  additional 
financial  support  through  grant  no.  0507-139. 

The  study  covered  by  this  report  was  under  the  direction  of  Professor 
E.  K.  Dabora ,  Department  of  Mechanical  Engineering  with  Mr.  James  J.  Murray 
of  ARO  as  contract  monitor  and  Dr.  Norman  Slagg  of  ARRADCCM  as  Scientific 
Liaison  Representative. 


ABSTRACT 


A  shock  tube  technique  was  used  to  determine  the  lean  detonation  limit 
of  sprays  of  kerosene  in  air  when  the  kerosene  is  sensitized  by  either 
propyl  nitrate  or  butyl  nitrite.  The  technique^r^described  in  detail  and -g — — 
is  used//as  a  first  step/^to  determine  the  relationship  between  the  detonation 
initiation  energy  and  the  initiation  power  for  stoichiometric  H?-air  mixtures 
at  .5  atmosphere.  It  is  found  that  if  the  initiation  energy,  E  is  expressed.. 


in  J/^n^  and  the  power  P,  in  Kw/jdn^  then  (E  -  10)  (P  -  27)  =  555 

method  is  also  used  to  check  on  the  lean  limit  of  H2-air  at 

f  atmospheric  pressure.  It  was  found  that,  for  ..an  initiation  energy  of  50  J/cm 

- ___  f  p  — 

at  P  =  160  Kw/£m^)  (initiation  M(=J3.5),  the  lean  limit  is  between  10-141  Hr 

in  air.  This  value  is  somewhat  lower  than  that  found  in  the  literature  whidry  — 


18. 


2%  H2*  in  air. 


Approximately  the  same  energy  and  power  were  used  to  determine  the  lean 

limit  ofL Jkerosene  sprays.  Monodisperse  sprays  having  droplets  diameter  = 
i'e  >  i  *  „  d 


780(jjp*were  used.  Three  equivalence  ratios  were  tested:  .59,  .  44/T.  3. 
Kerosene  was  mixed  with  10%  and  20%  of  either  propyl  nitrate  (PN)  or  butyl 
nitrite  (BN).  'Ther^results  indicate  that  the^addition  of  propyl  nitrate 

^  Cr^v>-  . 

reduces  the  limiting  equivalence  ratio Thus,  for  10%  PN  the  limiting  pj 


is  between  .44  and  .  3^and  for  20%  PN/  4;,  is  below  .3.  Erratic  behavior  was 
detected  when  10%  BN  was  used..  However,  when  20%  of  BN  was  mixed  in  kerosene 


the  limit ing/<t>\ was  again  below  .3. 
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I.  INTRODUCTION 


Two-phase  detonations  are  currently  of  interest  and  the  state  or 
their  understanding  has  been  reported  by  Dibora  and  Weinberger  (1)  and 
Nettleton  (2).  One  particular  type  of  two  phar.e  detonations,  namely 

i 

spray  detonations  has  been  investigated  in  some  depth  v/vthin  the  last 

I 

two  decades.  For  example,  the  work  of  Dabora  et  al.  (3)  and  Ragland  t 

et  al.  (4)  dealt  with  the  properties  and  structure  of  spray  detonations 
as  manifested  in  a  constant  area  tube  and  that  of  Nicholls  et  al.  (5) 
was  concerned  with  cylindrical  detonations.  A  more  up  to  date  work  in 
the  latter  category  is  given  by  Bar  Or  et  al.(6).  Lu  et  al.  (7)  were 
concerned  with  the  chemical  and  the  physical  processes  effects  on  the 
detonation  phenomenon  and  in  particular  found  that  certain  additives  to  the 
fuel  can  affect  the  establishment  of  detonations. 

More  recently  Dabora  (8,  9)  presented  a  theoretical  model  for  spray 
detonations  from  which  the  critical  energy  and  power  required  for  initiation 
could  be  established.  The  model  treats  all  three  types  of  detonations: 
planar,  cylindrical  and  spherical. 

Despite  the  rather  extensive  work  performed  so  far  on  spray  detonations, 
the  experimental  difficulties  involved  have  precluded  the  treatment  of  some 
aspects  of  the  phenomenon  in  an  exhaustive  manner.  For  exanple,  very  few 
systematic  studies  have  been  conducted  to  determine  the  detonability  limits. 
This  situation  is  to  be  contrasted  with  the  detonability  limits  of  gaseous 
mixtures  which  are  reasonably  well-known  for  many  different  combinations. 

Knowledge  of  detonability  limits  is  of  importance  from  at  least  two  con¬ 
siderations:  the  fuel-air  explosion  (FAE)  concept  and  hazards  elimination. 


In  the  former,  one  is  interested  in  conditions  that  would  insure  the 
occuranre  of  detonation  and  in  the  other, the  conditions  that,  would  reduce 
the  probability  of  detonations.  So  far  the  determination  of  limits  seems 
to  elude  theoretical  considerations  alone,  and  therefore,  almost  all  limits 
a.*e  based  on  either  direct  or  indirect  experiments. 

The  most  important  parameter  in  the  determination  of  limits  is  the 
initiation  energy.  For  detonable  mixtures  there  is  a  minimum  critical  energy 
that  is  usually  necessary  to  initiate  a  detonation,  and  for  a  given 
combination  of  fuel  and  oxidizer  the  minimum  critical  energy  depends  on  the 
mixture  ratio.  Near- stoichiome trie  mixtures  require  the  least  critical  energy 
but  as  the  mixture  becomes  either  fuel  rich  or  fuel  lean  the  critical  energy 
increases.  The  composition  limits  are  then  reasoned  to  be  ♦'hose  compositions 
that  require  inordinate  amounts  of  energy  before  a  detonation  would  occur. 

In  principle,  however,  true  limits  are  those  that  would  not  support  detonation 
no  matter  how  large  the  initiating  energy  is.  Limits  according  to  this 
definition  would  however  be  hard,  if  not  impossible  >to  determine,  a  more 
practical  approach  is  to  define  limits  as  those  compositions  that  require  an 
arbitrary  order  of  magnitude  higher  energy  than  say  the  stoichiometric  mixture 
would  require  for  a  detonation  to  take  place. 

Although  the  initiation  energy  seems  to  be  the  most  important  parameter, 
there  is  reason  to  believe  that  the  rate  of  energy  deposition  to  be  of 
importance  as  well  (10,  11).  Thus,  it  can  be  argued,  that  there  is  a  minimum 
initiation  energy  below  which  no  detonation  would  occur  no  matter  how  high 
the  power  is,  and  that  there  is  a  minimum  power  below  which  no  detonation  would 
occur  no  matter  how  high  the  total  energy  is.  This  idea  is  explored  somewhat 
in  this  report. 
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The  main  purpose  of  this  report  however,  is  to  determine  in  a  practical 
way  the  lean  limit  of  kerosene  spray  in  air,  when  two  additives  are  mixed 
with  kerosene:  namely,  propyl -nitrate  and  butyl -nitrate.  Since  one  would 
expect  that  the  limit  is  affected  by  drop  size., the  limit  must  be  determined 
for  each  drop  size.  Unfortunately  the  scope  of  the  work  was  limited 
because  of  time  consideration,  so  that  only  sprays  with  drop  size  of  the  order 
of  800  ym  were  considered. 
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1 1 .  EXPERIMENTAL  METHOD 


1.  Experimental  Apparatus: 

The  shock  tube  technique  was  used  to  determine  the  lean  limit.  In 
this  method  a  shock  wave  of  finite  strength  is  induced  into  the  driven 
section  of  the  shock  tube  which  is  prefilled  by  the  desired  mixture.  The 
decay  or  the  acceleration  of  the  wave  is  observed.  If  no  acceleration 
is  observed,  the  mixture  is  deemed  non-detonable.  If  however,  acceleration 
is  observed  which  is  indicative  of  a  detonation,  the  fuel  air  mixture  is  then 
leaned  in  subsequent  runs  until  no  detonation  takes  place  and  in  this  way 
the  lean  limit  is  identified. 

The  shock  tube  which  was  used  in  this  study  is  described  in  detail  in 
(12)  and  is  shown  schematically  in  Fig.  1  .  Briefly,  the  driven  section  of 
the  tube  is  9  ft.  long  with  a  2"  x  2"  square  cross-section.  The  tube  is 

topped  by  a  tee-section  to  accommodate  a  dump  tank  on  the  side.  Another  section 
on  top  of  the  tee  accommodates  a  spray  generator  which  is  similar  to  the 
one  described  by  Dabora  (13).  To  insure  a  reasonable  spray  along  the  tube, 
the  tube  is  oriented  in  the  vertical  direction.  The  driver  section  is  circular 
in  cross-section  providing  a  driver  to  driven  area  ratio,  Ai,/Ai  =  1.65.  The 
driver  length  is  variable  although  in  most  of  the  tests  a  6"  driver  length  was 
used.  The  tube  operates  on  the  double  diaphragm  principle  to  insure  diaphragm 
break  up  at  a  desired  pressure  ratio.  The  theoretically  obtainable  Mach 
numbers  are  shown  in  Fig.  2  when  the  driver  gas  is  helium  and  the  driven 
gas  is  air.  Generally  the  actual  Mach  numbers  obtained  were  closer  to  the 
Ah/Ai  ■  1  curve. 

The  tube  was  provided  with  pressure  switches  at  one  foot  intervals, 
starting  at  .5  ft.  from  the  diaphragm.  In  addition, two  piezoelectric  pressure 
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FIG.  1.  SCHEMATIC  OF  TUBE  AND  INSTRUMENTATION 
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transducers  are  placed  at  selected  locations.  At  location  6,  in  addition  to 
the  pressure  measurement,  optical  radiation  is  monitored  with  an  EGG  Lite 
Mike.  Photographic  observation  of  the  phenomenon  is  made  near  that  location 
with  a  spark  source  schlieren  system. 

The  wave  velocity  is  measured  from  time  intervals  as  recorded  in  a 
raster  display  of  the  pressure  switches  as  well  as  from  two  electronic  time 
interval  units.  Pressure  measurements  and  radiation  measurement  are  recorded 
with  oscilloscopes  as  indicated  in  Fig.  1. 

The  spray  runs  are  conducted  as  follows.  The  tube  is  purged  with  air, 
and  fresh  air  is  established  in  the  tube.  The  driver  and  the  diaphragm 
sections  are  filled  to  the  desired  pressure.  Through  a  mechanical  timer, 
the  drop  generator  is  started,  and  stopped  and  the  diaphragm  section  is  dumped 
to  the  atmosphere  via  a  solenoid  valve.  The  shock  is  thus  initiated  and  the 
measurements  are  automatically  recorded. 

2.  Theoretical  Background:  Shock  Wave  Decay: 

For  the  usual  use  of  a  shock  tube,  the  driver  length  is  chosen  long 
enough  to  insure  a  constant  wave  velocity  throughout  the  driven  section. 
However,  since  for  our  purpose  we  desired  to  observe  whether  a  detonation 
occurs  or  not  by  noting  the  wave  velocity,  it  was  decided  that  a  short 
driver  would  be  useful.  With  a  short  driver,  the  rarefaction  wave  in  the 
driver,  after  reflection  catches  up  with  the  wave  and  reduces  its  velocity 
if  the  wave  is  non- reactive.  If  the  wave  is  reactive  as  in  a  detonation,  it 
i.i  expected  that  the  wave  will  continue  to  propagate  unhampered  by  the 
rarefaction  wave. 


In  order  to  get  an  approximate  value  of  the  decay  of  the  wave,  we  note 

that  the  initiation  Mach  number  of  interest  is  between  3  and  4.  For  the  gases 

used,  helium  in  the  driver  and  air  in  the  driven  section,  the  Mach  number 

obtainable  for  a  given  pressure  ratio  is  shown  in  Fig.  2.  Fox  (12)  has 

calculated  the  location  at  which  the  rarefaction  wave  catches  up  with 

shock  wave  as  a  function  of  the  driver  length.  The  ratio  of  this  critical 

location  over  the  driver  length,  X  /f4,  is  found  to  be  2. 8-5. 2  for 

>»• 

M  =  3-4  when  the  area  ratio  is  1.65.  If  after  that  location  the  wave  is 
considered  as  a  constant  energy  blast  wave  then  its  decay  can  be  approximated 


(14)  by 

dM  a+1  n  1  s  dR 

H - r  V'W)  TT 


(1) 


where  R  is  the  blast  wave  distance  and  a  is  a  geometric  factor  =  0,  1,  2 
for  planar  .cylindrical  or  spherical  wave.  For  the  shock  tube,  a  =  0  and 
therefore: 

p) 


This  equation  can  be  integrated  to  give: 

R  _  Mo2-1 
Ko  "  FTT“ 


(3) 


For  a  20%  decay  in  M,  R/R0  =  1.65  for  M0  ■  3-5  and  for  a  50%  decay,  R/R0 
varies  from  4.57  to  6.40  when  M0  varies  from  5  to  3.  If  R©  is  assumed 
to  be  the  location  where  the  rarefaction  overtakes  the  shock  wave,  the  driver 
length  can  then  be  designed  so  that  observable  decay  could  be  effected 
within  a  reasonable  length  of  the  driven  section.  Thus  for  the  6"  driver 
at  Mo  ■  3.5,  for  which  X  /-€.«♦  ■  Ro/£it  ■  4,  20%  decrement  is  expected  at 
.5  x  4  x  1.65  ■  3.3  ft.  from  the  diaphragm. 
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Some  tests  were  conducted  to  check  the  predicted  velocity  decrement. 

For  a  6"  long  driver  the  velocity  decrement  of  a  shock  wave  in  air  at 
M  2  3.6  resulted  in  a  decrement  of  17%  at  6  ft.  from  the  diaphragm.  Thus 
it  appears  that  the  experimental  decrement  is  about  half  the  theoretically 
predictable  decrement. 

3,  The  Relation  Between  Energy  and  Power  for  Detonation  Initiation: 

It  was  indicated  in  the  introduction  that  the  relationship  between  power 
and  energy  for  detonation  initiations  would  be  examined.  In  a  shock  tube 
experiment  the  power  is  effectively  related  to  the  initiation  Mach  number  and 
the  energy  is  related  to  the  effective  energy  stored  in  the  driver.  Thus 
for  a  given  set  of  driver  and  driven  gases,  the  power  is  related  to  the  pressure 
ratio  pu/pi  and  the  energy  to  the  length  of  the  driver. 

The  rate  of  energy  input  to  the  driven  section  must  equal  the  rate  of 
work  performed  by  the  interface.  Thus  the  rate  per  unit  area  is  =  p2U2,  or 
the  power  per  unit  area  is: 

P  -  P2U2  (4) 

=  Pi(P2/pO(ui-u2) 

*  (P2/P1HI  -  U2/Ui)piUi 

*  (P2/P1HI  -  Pi/p2)piMiai  (5) 

Eq.  5  can  be  written  in  terms  of  the  shock  Mach  numbers  after  using  the 
shock  relations.  Thus: 

P  -  * 1]  tfffnnft  P’ai  (6) 

The  energy  involved  is  equal  to  this  rate  multiplied  by  the  time  in  which 
this  rate  is  effective.  Hus  time  is  considered  to  be  the  time  in  which  the 
wave  propagates  until  the  rarefaction  wave  reflected  from  the  driver  end  wall 
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catches -up  with  the  shock  front.  Thus  if  Xc  is  the  distance  from  the 
diaphragm  at  which  the  rarefaction  wave  reaches  the  shock  front,  then 
the  energy  per  unit  area: 


£  *  P  Xc/u, 

(7) 

or 

E  =  P(Xc/£4K4/Miai 

(8) 

where  £4  is  the  driver  length  and  Xc/£4  depends  on  Mi  for  a  given 
combination  of  gases. 

Some  tests  on  the  initiation  of  Ife-air  mixtures  were  performed. 

As  background,  the  detonation  velocity  and  Mach  number  have  been  calculated 
using  the  computer  code  of  Gordon  and  McBride  (15).  These  are  shown  in 
Fig.  3  and  4.  The  initiation  of  a  stoichiometric  mixture  initially  at  a 
pressure  of  0.5  atmosphere  was  checked  with  three  drivers  of  different 
lengths  (2",  4"  and  8M)  with  varied  pressure  ratios.  The  location  at  which  the 
Chapman- Jouguet  (CJ)  Mach  number  is  reached  is  noted  and  the  results  are 
shown  in  Fig.  5.  It  can  be  seen  that  the  8"  and  4"  results  show  reasonable 
variations  whereas  the  2"  results  show  much  scatter.  At  any  rate,  the  pressure 
ratio  below  which  detonation  is  expected  to  take  place  at  extremely  long 
distance  from  the  diaphragm  can  be  estimated.  The  corresponding  initiation 
Mach  number  can  be  found  from  Fig.  2  and  the  corresponding  Xc/£ 4  is  found 
from  (12).  Equations  (6)  and  (8)  are  then  used  to  determine  the  power  and 
energy  for  initiation.  The  results  are  presented  in  Table  I ,  and  plotted  in 
Fig.  6.  From  this  figure  the  asymptotes  of  energy  and  power  can  be  estimated. 
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Table  I 


ENERGY  AND  POWER  RELATION  FOR  THE  INITIATION  OF 
STOICHIOMETRIC  Hj-AIR  AT  -  0.5  atm. 


Driver  Length 
Ph (in) 

p*/pi 

Mi 

XC/*H 

k 

P 

watts/cm2 

E 

J/cm2 

2 

100 

3.90 

4.8 

24.4 

10.89  x  10* 

16.69 

4 

65 

3.55 

3.8 

38.6 

8.08  x  10* 

21.50 

8 

40 

3.2 

3.1 

63 

5.78  x  10* 

27.89 

These  would  correspond  to  the  minimum  initiation  energy  and  power,  ihe 
minimun  energy  can  also  be  determined  by  plotting  E  vs.  Xc,  as  is  done  in 
Fig.  7,  and  finding  the  internet.  This  gives  E.n  ■  10  J/cm2  for  a 
stoichiometric  mixture  of  H2-air  at  0.5  atmos.  If  the  induction  distance  is 
is  inversely  proportional  to  pressure  (16) ,  the  minimun  initiation  energy 
at  atmospheric  pressure  becomes  5  J/cm2  since  the  energy  for  a  planar 
wave  can  be  considered  proportional  to  the  induction  distance. 

It  can  be  shown  from  a  plot  of  F  vs.  1/X„  that  P  •_  «  2.7  x  10“  watts/cm2 

c  min 

which  corresponds  to  an  initiation  Mach  nunber  of  2.55  The  two  asynptjtic 
values  of  and  P^  can  be  combined  to  give  the  following  relation 
between  P  and  E: 

(P  -  2.7  x  10*)  (E  -  10)  -  55.5  x  10*  (9) 

This  equation  is  at  variance  with  the  relationship  of  energy  and  power 
found  by  Khystautas  and  Lee  (11)  for  spark  initiation  of  stoichiometric 
acetylene- oxygen  detonations.  In  their  case  they  found  a  double  value  of 
energy  for  a  given  power  which  appears  to  be  physically  unsatisfactory. 

For  the  question  arises  as  to  why,  for  a  given  initiating  source  operating 
at  a  certain  rate  of  energy  deposition,  there  should  result  a  detonation  at 
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two  values  of  energy  deposition,  one  high,  one  low.  If  the  energy  deposition 

could  be  stopped  at  the  lower  value,  a  detonation  should  still  be  possible, 

be  initiated  and  the  additional  input  energy  would  seem  superfluous.  It  is 

hoped  that  future  work  would  resolve  this  apparent  anomaly. 

At  any  rate,  the  values  of  R  .  and  P.„  found  here  can  be  used  as 

nun  min 

guidelines  for  the  determination  of  limits.  The  lean  limit  of  Hj-air  will 
be  checked  first. 
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III.  THE  LEAN  LIMIT  OF  H»-AIR 

It  was  thought  that  the  method  for  determining  the  limits  of 
kerosene-air  should  first  be  tried  for  a  combination  of  fuel-oxidizer  of 
known  limits.  Since  the  detonation  limits  of  Ha -air  are  supposedly  well 
known,  this  mixture  was  chosen  for  testing.  The  accepted  value  of  the 
lean  limit  is  18.21  Ha  in  air  at  0  *  .S3  (17,  18). 

Tests  were  made  on  Ha -air  mixtures  at  0,  10,  14,  16,  20  and  301  Ha 
to  determine  the  initiation  behavior.  The  mixtures  were  at  an  initial 
pressure  of  1  atmosphere  and  the  initiation  Mach  number  was  about  3.5. 

The  driver  length  was  6"  providing  an  approximate  energy  input  50  J/cm2 
at  a  power  of  16  x  10**  watts/cm2.  The  energy  input  is  about  an  order  of 
magnitude  larger  than  what  is  needed  to  initiate  the  stoichiometric  mixture. 

For  each  run, the  velocity  was  monitored  throughout  the  tube  and  pressure 
measurements  were  made.  In  particular  the  pressure  at  location  6  (Fig.  1), 

5  1/2  ft.  fromrthe  diaphragm  as  well  as  light  emission  were  recorded. 

Figure  8  shows  representative  oscilloscope  records  for  the  pure  air  case, 
the  101  H2  and  the  201  H2  runs.  The  pressure  records  seem  to  be  identical 
in  shape.  This  reflects  the  fact  that  whether  an  adiabatic  shock  wave  or 
a  detonation  wave  is  involved,  the  wave  is  followed  by  an  expansion.  In  the 
detonation  case  this  is  an  inherent  phenomenon.  In  the  shock  case,  this  is 
due  to  the  limited  driver  length  and  the  resultant  early  interaction  of  the 
reflected  expansion  wave  with  the  shock  wave. 

The  lack  of  emission  in  the  pure  air  is, of  course,  expected.  The  lack 
of  emission  for  the  10%  Ha  is  taken  as  indication  of  the  lack  of  reaction  and 
therefore,  absence  of  detonation.  In  the  201  H2  case,  emission  is  prominent  and 
indicative  of  at  least  some  reaction. 


18 


Run  #54 


Trig,  delay  from  3  =  dO 
llorz.  =  100  ps/div 
Vert.  Top  =  10  mV/div 

Hot  =  100  psi/div 


(h)  10°o  11; 

Run  #03 

Trig,  delay  from  3  =  50 
llorz.  =  200  us/d iv 
Vert.  Top  =  5  mV/div 

Rot  =  100  psi/div 


(cj  20*  Hg 
Run  *49 

Trig,  delay  from  3  =  25 
llorz.  =  100  ps/di v 
Vert.  Top  =  20  mV/div 

Rot  -  200  psi/div 


no.  s  i:missio\  ,v\d  prhssurj:  Rirouns  \t  rotation  o  tor  air,  air 

10':,  !i:  and  AIR  +  201  Hr 


Figure  9  shows  the  Mach  nunber  variation  throughout  the  tube  for  the 
six  cases.  It  is  seen  that  for  the  pure  air  case  and  the  10%  II2  case  the 
wave  slows  down  generally.  In  the  14%  H2  case  the  Mach  number  can  reach 
a  value  of  3.9  which  is  the  corresponding  CJ  Mach  number.  In  the  16%, 

20%  and  30%  H2  cases  the  Mach  number  exceeds  the  corresponding  CJ  Mach 
number  of  4.12,  4.45  and  4.85  respectively.  However,  no  steady  detonation 
could  be  observed  as  the  Mach  nunber  seems  to  decrease  after  a  peak.  It 
is  felt  that  if  the  tube  length  could  be  increased  a  steady  detonation  would 
eventu&lly  develop.  Unfortunately  it  was  not  possible  to  do  this.  However, 
some  basis  for  determining  the  limit  had  to  be  made. 

Ttoo  criteria  were  used;  one,  indication  of  emission  and  the  other, the 
attainment  at  same  point  in  the  tube  of  a  Mach  nunber  squal  or  exceeding  the 
CJ  Mach,  especially  when  the  initiating  Mach  nunber  is  lower  than  the  CJ 
Mach  nunber  as  in  these  runs.  Based  on  these  two  criteria,  it  appears  that 
the  lean  detonation  limit  for  H2-air  as  determined  by  our  method  is  somewhere 
between  10%  H2  and  14%  H2  which  is  loweT  than  the  currently  accepted  value. 

At  this  point,  it  is  appropriate  to  mention  that  recent  experiments (19) 
on  shock  initiation  of  H2-02  mixtures  with  argon  dilution  indicate  that 
conditions  exist  in  which,  instead  of  obtaining  CJ  detonations,  "reaction 
enhanced"  waves  are  observed.  These  waves  are  observed  in  weakly  reactive 
mixtures  (due  to  either  low  reactants'  concentrations  or  high  inert  dilution) 
wherein  only  part  of  the  possible  heat  release  influences  the  wave  velocity. 
The  experiments  were  conducted  with  a  long  driver  such  that  the  adiabatic 
wave  would  propagate  without  decay.  In  contrast,  our  experiments  are  at  high 
reactants'  concentrations  (high  pressure)  and  in  a  tube  with  a  short  driver. 
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FIG,  9,  PROPAGATION  MACH  NUPBER  VARIATION  WITH  DISTANCE 


In  the  experiments  of  (19)  "the  reaction  enhanced"  waves  do  not  reach 
or  exceed  the  CJ  velocity  when  the  initiation  Mach  number  is  below  the  CJ 
Mach  number.  However,  in  our  experiments  we  detect  regions  in  which  the 
CJ  velocity  is  exceeded.  Thus ,  it  is  felt  that  our  method  which  is  based 
on  a  decaying  adiabatic  wave  is  a  reasonable  method  for  determining  limits. 


22 


IV.  KEROSENE-AIR  EXPERIMENTS 


1.  Theoretical  Calculations: 

The  Chapman  Jouguet  detonation  properties  of  mixtures  of  kerosene- air 
and  kerosene  +  additives-air  were  calculated  using  the  computer  code  of 
Gordon  and  McBride  (15).  Kerosene  was  assumed  to  have  an  average  formula 
of  Ci 1.6  H2S.2  with  a  heat  of  formation  of  -5.6  Kcal/mole.  The  additives 
used  were  n-propyl  nitrate  (C3H7NO3)  and  n-butyl  nitrate  (C4H9NO2).  Their 
heats  of  formation  were  taken  as  -41.6  and  -50.0  Kcal/mole.  In  the  crlculations 
the  fuel  was  assumed  to  be  in  the  gaseous  phase,  an  assumption  which  has 
practically  no  effect  on  the  detonation  velocity.  However,  the  fact  that 
the  fuel  is  liquid  does  have  an  effect  on  the  detonation  Mach  number  and 
the  pressure  ratio.  The  effect  depends  on  the  loading  factor  or  the  mass 
ratio  of  the  liquid  component  (20).  For  lean  mixtures  which  are  of  interest 
in  this  work  the  effect  is  negligible. 

The  calculated  properties,  such  as  the  detonation  velocity,  Mach  nunber 
and  pressure  ratio  are  plotted  in  Figs.  (10,  11,  12)  respectively.  These 
plots  do  not  include  any  correction  due  to  the  loading  factor. 

2.  Fuel-Air  Mixture  Ratio  Determination: 

The  tests  were  to  be  made  in  monodisperse  sprays.  As  was  indicated 
before  the  spray  is  formed  by  a  method  similar  to  that  of  Dabora  (13)  where 
capillaiy  jets  are  broken  up  into  equal  size  droplets  by  inducing  deliberate 
disturbances  into  the  liquid  jets  at  an  appropriate  frequency.  The  nominal 
size  of  the  drops  desired  was  700  pm.  The  closest  capillary  needle  that  can 
produce  this  drop  size  is  guage  #22. with  an  internal  diameter  of  .413  nm.  The 
drop  size  produced  is  .413  x  1.89  ■  .780  nm  or  780  pm. 


TION  VELOCITY  -  M/SEC 


1900 


0  .50  1.0  1.5  2.0  2.5 

EQUIVALENCE  RATIO 

FIG,  10,  CHAPMAN  JOJGUET  DETONATION  VELOCITY  OF  KEROSENE, 
WITHOUT  AND  WITH  ADDITIVES,  AND  AIR 
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0  .5  1.0  1.5  2.0  2.5  3.0 

EQUIVALENCE  RATIO 

FIG,  11,  DETONATION  MACH  NLfBER  OF  KEROSENE  (+  ADDITIVES)  AND  AIR 
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EQUIVALENCE  RATIO 

FIG.  12,  DETONATION  PRESSURE  RATIO  OF  KEROSENE  (+  ADDITIVES)  AND  AIR 
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According  to  (13)  the  appropriate  frequency  and  the  shedding  velocity 
of  a  kerosene  jet  are  390  Hz  and  74  an/sec.  The  terminal  velocity  of  the 


droplet  is  250  an/sec. 

The  air- fuel  ratio  per  single  colunn  of  fuel  drops  can  be  calculated 
from  the  following  equation: 


m 


air 


m 


fuel 


6  *Valr 

TTd*  f  p £ 


(10) 


where  A  is  the  area  of  the  shock  tube,  u^  the  terminal  velocity  of  the 
droplets,  d,  their  diameter  and  f,  the  frequency.  Using  the  appropriate 
values  (A  -  25.8  cm2),  the  air  fuel  ratio  is  found  to  be  99.6  per  colunn  of 
drops.  The  stoichiometric  air- fuel  ratio  as  determined  from  the  stoichiometric 
equation  is  found  to  be  14.76.  So  that  the  equivalence  ratio  per  column  is 
0.148.  Because  of  the  manner  in  which  the  spray  is  generated  the  equivalence  ratio 
can  be  varied  in  a  discrete  manner  only.  Thus  equivalence  ratios  of  .3,  .44, 

.59,  etc.  are  obtainable. 

It  should  be  mentioned  that  oecause  the  shedding  velocity  is  smaller  than 
the  terminal  velocity, the  fuel-air  ratio  is  larger  near  the  top  of  the  tube 
than  the  bottom  of  the  tube.  It  is  possible  to  estimate  the  distance  for 
the  drops  to  travel  before  they  reach  a  velocity  within  951  of  the  terminal 
velocity  (20).  For  our  conditions  this  turns  out  to  be  about  3  ft. 

In  the  calculations  above  the  effect  of  propyl -nitrate  and  butyl-nitrite 
is  assumed  to  be  negligible.  The  reason  is  that  the  maximun  amounts  used  are 
201  by  weight  and  their  densities  are  close  to  that  of  kerosene.  Furthermore, 
the  carbon  to  hydrogen  ratio  is  approximately  the  same  as  that  of  kerosene. 


3.  Energy  Requirement  for  Initiation: 

The  minimum  energy  requirement  for  the  initiation  of  detonations  in 
sprays  depends  on  many  factors  such  as  the  drop  size,  the  mixture  ratio,  the 
ignition  delay  and  the  energy  deposition  profile  after  ignition  (9).  The 
last  two  properties  are  not  known  for  either  kerosene  or  kerosene  with 
additives,  and  it  is  therefore,  necessary  to  make  some  assumptions  before  the 
minimum  initiation  energy  could  be  estimated.  If  the  ignition  delay  is 
assumed  to  be  equal  to  one-half  the  drop  breakup  time  and  if  the  heat  release 
from  each  drop  is  assumed  to  be  instantaneous,  then  the  method  of  (9)  can 
provide  an  estimate  of  the  minimum  initiation  energy.  From  Table  1  and 
Fig.  2  of  (9) ,  it  is  found  that  E*/d  =  100  J/cm*  which,  for  d  *  780  u,  gives 
E*  *  7.8  J/cm2.  As  indicated  in  Section  III,  a  6"  driver  provides  “50  J/cma 
when  the  initial  Mach  number  is  3.5.  It  is  also  indicated  in  (9)  that  if  the 
breakup  time  is  controlling  then  the  minimum  initiation  energy  is  minimally 
affected  by  the  equivalence  ratio.  On  the  other  hand,  if  the  ignition  delay 
is  controlled  by  a  chemical  induction  time  the  minimum  ignition  energy  can  be 
two  orders  uf  magnitude  larger  for  4>  =  .5  than  that  of  <j»  *  1.  Thus  the 
50  J/cm2  provided  by  the  driver  used  in  our  experiments  can  be  considered 
adequate  for  the  determination  of  limits, only  if  the  drop  breakup  time  is 
controlling  the  ignition  delay. 

4.  Experimental  Results: 

Previous  work  (21)  has  indicated  that  kerosene-air  mixtures  are  non- 
detanable.  Limited  runs  conducted  here  with  pure  kerosene  at  $  *  .44  produced 
no  detonation  either.  Thus,  our  work  concentrated  on  the  effect  of  the  two 
additives  -  propyl  nitrate  and  butyl  nitrite. 
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Specifically  tests  at  $  *  .59,  .44  and  .3  were  conducted  for  101  and 
201  of  the  additives  in  the  fuel.  For  each  test  a  schlieren  photograph 
was  taken  to  insure  that  the  spray  is  present  and  to  give  us  a  visual  feel 
of  the  phenomenon.  In  addition,  pressure  data  and  emission  records  were 
taken.  Also,  the  progress  of  the  wave  is  monitored  throughout  the  tube.  Not 
all  of  the  instrunentatian  worked  perfectly  every  run,  so  that  some  runs 
had  to  be  discarded  because  of  incomplete  data.  As  indicated  before,  the 
criteria  for  detonation  were  light  emission  and  the  attainment  of  a  wave 
speed  equal  or  higher  than  the  CJ  wave  speed  at  some  point  in  the  detonation 
tube. 

Representative  data  are  reported  here.  For  the  101  propyl  nitrate, 

Figs.  (13  §  14)  show  schlieren  photographs,  pressure  and  emission  records 
at  <J>  =  .44  and  .3.  The  schlieren  photograph  in  Fig.  13  shows  a  typical 
detonation  wave  in  which  the  front  appears  as  a  non-planar  front  because  of 
some  interactions  from  waves  behind  it.  The  pressure  record  shows  some 
oscillations  and  pronounced  peaks  which  could  be  due  to  blast  waves  from 
each  drop  although  such  blast  waves  do  not  seen  to  be  detectable  in  the  schlieren 
photograph.  The  emission  record  seems  to  indicate  the  presence  of  the  peaks 
also.  One  interesting  feature  is  the  very  short  ignition  delay. 

The  schlieren  photograph  in  Fig.  14  shows  a  typical  non-reactive  front. 
Usually  in  such  a  situation  the  front  remains  planar.  The  pressure  record  shows 
no  pronounced  peaks  and  the  emission  trace  shows  no  radiation. 

The  wave  velocity  variation  is  shown  in  Fig.  15.  It  can  be  seen  that  for 
4>  «  .59  and  .44  the  propagation  velocity  increases  at  some  point  in  the  tube 
and  reaches  a  velocity  conparable  to  the  respective  CJ  velocity.  On  the  other 
hand  for  $  ■  .3  the  velocity  does  not  veer  too  much  from  the  initiation  velocity 
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Run  #130 


Trig,  delay  from  3  =  400  nsec 
Ilorz.  =  200  ns/di v 
Vert.  Top  =100  psi/div 
Rot  =  5  mV/di v 


FIG.  13  SGIUEEEN  PI IOTOGRAPI I ,  EMISSION  AND  PRESSURE  RECORD 
FOR  KEROSENE  +  Utt  PROPYL  NITRATE  AT  LOCATION  6  (<f>  =  .44) 
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Run  #128 

Tube  width  =  5.08  cm* 


Trig,  delay  from  3  =  400  us 
Morz.  =  200  us/div 
Vert.  Top  =  100  psi/div 
Bot  =  5  mV/div 


FIG.  14  SQ1LIEREN  PHOTOGRAPH,  EMISSION  A.‘JD  PRESSURE  RECORD 
FOR  KEROSENE  +  10”,  PROPYL  NITRATE  AT  LOCATION  6  (<p  =  .3) 
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DISTANCE  FROM  DIAPHRAGM  -  FT 

FIG,  15  VARIATION  OF  PROPAGATION  MACH  NUMBER  WITH  DISTANCE  FOR  KEROSENE  +  10%  PROPYL-NITRATE 


and  eventually  drops  down  to  below  the  CJ  velocity.  This  situation  together 
with  the  emission  records  is  interpreted  as  absence  of  a  detonation.  Thus, 
the  conclusion  is  that  the  lean  limit  for  kerosene  with  101  propyl -nitrate 
is  somewhere  between  $  *  .44  and  .3.  As  indicated  before,  because  of  our 
method  of  generating  the  spray,  it  is  not  possible  to  give  a  more  precise 
limit. 

The  results  on  runs  with  201  propyl  nitrate  are  exemplified  by  Figs. 

16  and  17.  For  the  (j>  *  ,44  case,  Fig.  16,  similai  conments  to  the  10% 
propyl  nitrate  case  apply.  For  the  $  *  .3  case,  Fig.  17,  the  wave  again  is 
non  planar  and  the  pressure  record  shows  peaks.  Hie  emission  trace  shows 
that  a  reaction  has  taken  place,  however  the  rise  in  emission  takes  longer 
time  than  for  the  <J>  *  .44.  The  velocity  data  of  Fig.  18  show  that  a  detonation 
is  possible  down  to  <}>  ■  ,3  and  therefore,  the  detonation  limit  in  this 
case  is  lower  than  <J>  «  .3.  Thus  it  appears  that  propyl  nitrate  does  affect 
the  limit  in  a  qualitative  way:  the  higher  its  mass  fraction  in  the  fuel, 
the  lower  the  lean  limit. 

The  results  with  101  butyl  nitrate  were  rather  erratic.  There  seemed  to 
be  cases  of  no  detonation  even  at  $  ■  .59  whereas  sometimes  detonations  at 
$  *  .3  could  be  observed.  As  a  result  of  this  situation  the  limit  for  this 
mixture  could  not  be  determined  with  any  certainty. 

For  the  201  butyl  nitrite  mixture,  the  data  seem  to  be  comparable  to  the 
20%  propyl  nitrate.  Figure  19  shows  typical  results  for  <t>  -  .44  and  Fig.  20 
is  indicative  of  the  results  for  d>  ■  .3.  Again  detonation  appears  possible 
at  <J>  ■  .3  as  cirribirated  by  the  wave  velocity  data  shown  in  Fig.  21.  Thus 
the  limit  for  this  mixture  is  below  $  ■  .3. 
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Run  #121 


<? 


Tube  width  =  5.08  cm. 


Trig,  delay  from  3  =  400  ysec 
Horz.  =100  ys/div 
Vert.  Top  =  200  psi/div 
Bot  =  10  mV/div 


FIG.  16  SCHLIEREN  PHOTOGRAPH,  EMISSION  .AND  PRESSURE  RECORD 
FOR  KEROSENE  +  20!  PROPYL  NITRATE  AT  LOCATION  6  (4>  =  .44) 


34 


Run  #121 

Tube  width  =  5.08  cm. 


Trig,  delay  from  3  =  400  usee 
Iforz.  =  100  us/div 
Vert.  Top  =  200  psi/div 
Rot  =10  mV/div 


PI(I.  10  SUIMEREN  PIIOTiXIRM’H,  EMISSION  .AND  PRESSURE  RECORD 

i :or  kerosene  +  2o°0  propyl  nitrate  at  location  o  (<j>  =  .44) 
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Tube  width  =  5.08  cm. 


(<•  17  sciiueren  piiotoorapii,  emission  and  pri-ssurt  record 
I'OR  KEROSENE  +  20  o  P ROPY  I,  NITRATE  AT  LOCATION  6  (d>  -  .3) 


DISTANCE  FROM  DIAPHRAGM  -  FT 

VARIATION  OF  PROPAGATION  MACH  NUMBER  WITH  DISTANCE  FOR  KEROSENE  +  20%  PROPYL-NITRATE 


Run  #145 


Tube  width  =  5.08  cm. 


Trig,  delay  from  3  -  401 
Horz.  =  200  ps/div 
Vert.  Top  =  100  psi/div 
Bot  =  5  mV/div 


PRESSURE  RECORD 


6  =  .44) 


Run  #148 


Trij>. 

delay 

from  7>  =  400 

Ilorz. 

=  200 

ps/di V 

Vo  rt . 

Top  = 

100  psi/div 

Hot  = 

S  mV/div 

10.  20  SQILIEREN  PHOTOGRAPH,  OMISSION  AND  Plil '.SSIIRH  RECORD 
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In  brief  then,  at  an  initiation  energy  of  approximately  50  J/an*  the 
lean  limit  of  kerosene  with  101  propyl  nitrate  in  air  is  between  $  •  .3  and 

♦  *  .44.  For  201  of  propyl  nitrate  or  butyl  nitrate  the  limit  is  below 

♦  *  .3. 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  shock  initiation  technique  is  shown  to  be  an  adequate  technique 
to  determine  the  limits  for  direct  initiation  of  detonation  waves.  The 
technique  rests  on  the  idea  that  if  no  reaction  takes  place,  the  wave  would 
slow  down,  provided  a  proper  length  driver  is  chosen.  On  the  other  hand,  if 
a  detonation  occurs  the  wave  velocity  would  increase  or  at  least  the  wave 
will  be  self-driven  after  the  effect  of  the  driver  ceases. 

The  shock  initiation  technique  can  be  used  to  determine  the  relationship 
between  energy  and  power  necessary  to  induce  direct  detonation.  This  was 
done  for  stoichiometric  H2-Air  mixtures  and  it  was  found  that  in  addition  to  a 
miniimir.  energy  requirement,  there  is  a  minimum  power  requirement  as  well. 

For  values  higher  than  these  minima  a  monotonic  relationship  between  power 
and  energy  exists. 

The  technique  was  used  to  check  the  lean  limit  of  H2-Air.  The  limiting 
equivalence  ratio  was  found  to  be  lower  than  that  quoted  in  the  literature. 
This  suggests , as  per  current  awareness  in  the  literature,  that  the  method  of 
initiation  can  have  a  definite  effect  on  the  initiation  energy  and  therefore, 
on  the  limiting  equivalence  ratio.  The  technique  needs  to  be  explored 
further  for  other  combustible  mixtures,  and  the  results  need  to  be  compared 
with  those  of  other  initiation  methods  ascertain  the  reasons  for  any 
differences. 

Although  the  shock  tube  technique  was  used  here  for  spray  detonations, 
it  should  be  realized  that  some  fundamental  data  on  the  ignition  of  sprays 
is  still  necessary  before  the  technique  could  be  considered  completely  valid. 
In  particular  ,data  on  ignition  delays  and  comparison  with  drop  breakup  delays 
are  needed.  The  ignition  delay  affects  the  miniimjn  initiation  energy  which  in 
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turn  can  be  used  as  a  guideline  for  what  energy  input  should  be  used  to 
detenrdne  the  conposition  limit. 

The  composition  limit  for  kerosene  with  propyl  nitrate  and  butyl  nitrite 
was  determined  here.  However,  the  energy  input  was  limited  to  50  J/cm*. 
and  the  drop  size  to  800  urn.  Other  drop  sizes  need  to  be  investigated  and 
the  effect  of  larger  energy  input  should  be  ascertained. 

The  technique  used  here  is  limited  to  planar  geometry.  The  application 
of  the  results  to  other  geometries  in  particular  to  the  spherical  geometry 
(mean fined  detonations),  need  to  be  checked.  Pmdamental  data  such  as 
ignition  delay  and  energy  deposition  rates  are  deemed  necessary  to  carry  over 
the  planar  results  to  the  spherical  geometry.  A  relatively  inexpensive  manner 
to  effect  this  is  to  use  gaseous  detonative  mixtures  which  offer  handling 
simplicity. 
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VI.  NOMENCLATURE 


Ai  -  Driven  section  area 

Aj  »  Driver  section  area 

ai  *  Speed  of  sound  of  gas  in  driven  section 

d  -  Drop  diameter 

E  *  Energy 

E*  -  Critical  initiation  energy 
f  -  Droplet  shedding  frequency 
Is  *  Driver  length 
M  *  Mach  number 
M0  *  Initial  Mach  number 
m  *  Mass 

P  *  Power 

pi  *  Original  pressure  in  driven  section 

pa  *  Pressure  behind  shock  wave 

Ph  *  Driver  pressure 

R  »  Distance  from  initiation  source 

R  *  Gas  constant 

Ui  *  Shock  velocity 

u2  *  Velocity  downstream  of  shock 

ua  ■  Convective  velocity  behind  shock  wave 

U£  *  Terminal  velocity  of  droplet 

X  *  Distance  from  diaphragm  at  which  the  rarefaction  wave 
c  interacts  with  shock  front. 

a  ■  Geometric  factor  -  0,  1,  2  for  planar,  cylindrical  or 
spherical  wave  respectively 

*  Equivalence  ratio  ■  fuel-air  ratio/ (fuel -air  ratio)  stoich. 


p  «  Density 

Px  ■  Density  ahead  of  wave 
Pa  ■  Density  behind  wave 
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A  shock  tifce  technique  was  used  to  determine  the  lean  detonation  lumt 
,  nr  kgrosime  in  air  when  the  kerosene  is  sensitized  by  either  propyl 

nitrate  or  butyl  nitrite.  The  technique  is  described  ^  used» 

«  «  ««♦  van  to  determine  the  relationship  between  the  detonation 
initiation  energy  and  the  initiation  power  for  stoichiometric  H2-air  mixtures 
2*  i  ISLSSm  It  isfAnd  that  if  the -initiation  energy,  E  is  expressed 
2  i/cTS^'po^r then  (E  -  10)(P  -  27)  -  555  J-KW/cm>. 
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Abstract  Continued. 


the  aethod  is  also  used  to  check  on  the  lean  limit  of  Hi -air  at 
atmospheric  pressure.  It  was  found  that  for  an  initiation  energy  of  SO  J/ca2 
at  P  *  160  KW/ca2  (initiation  M  -  3.5),  the  lean  limit  is  between  10*14%  H2 
in  air.  This  value  is  somewhat  lower  than  that  fotnd  in  the  literature  which 
is  18.21  Hi  in  air. 

Approximately  the  same  energy  and  power  were  used  to  determine  the  lean 
limit  of  kerosene  sprays.  Monodisperse  sprays  having  droplets  diameter  * 

780  im  were  used.  Three  equivalence  ratios  were  tested:  .59,  .44,  .3. 

Kerosene  was  mixed  with  10%  and  20%  of  either  propyl  nitrate  (PN)  or  butyl 
nitrite  (BN).  The  results  indicate  that  the  addition  of  propyl  nitrate  reduces 
the  limiting  equivalence  ratio,  4.  Thus,  for  10%  PN  the  limiting  4  is  between 
.44  and  .3  and  for  20%  PN,  4  is  below  .3.  Erratic  behavior  was  detected  when 
10%  BN  was  used.  However,  when  20%  of  BN  was  nixed  in  kerosene  the  limiting 
4  was  again  below  . 3. 


